Continuous-time quantum walks on Erdos-Renyi networks 
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We study the coherent exciton transport of continuous-time quantum walks (CTQWs) on Erdos- 
Renyi networks. The Erdos-Renyi network of A*' nodes is constructed by connecting every pair of 
nodes with probabihty p. We numerically calculate the ensemble averaged transition probability of 
quantum transport between two nodes of the networks. For finite networks, we find that the limiting 
transition probability is reached very quickly. For infinite networks whose spectral density follows 
the semicircle law, the efficiencies of the classical and quantum-mechanical transport are compared 
on networks of different average degree k. In the long time limiting, we consider the distribution 
of the ensemble averaged transition probabilities, and show that there is a high probability to 
find the exciton at the initial node. Such high return probability almost do not alter in a wide 
range of connection probability p but increases rapidly when the network approaches to be fully 
connected. For networks whose topology is not extremely connected, the return probability is 
inversely proportional to the network size A'^. Furthermore, the transport dynamics are compared 
with that on a random graph model in which the degree of each node equals to the average degree 
k of the Erdos-Renyi networks. 

PACS numbers: 05.60.Gg, 03.67.-a, 05.40.-a 
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I. INTRODUCTION 

During the last few years, the coherent exciton dynam- 
ics in quantum system has been extensively studied hy 
both experimental and theoretical methods [11, H, H, [J. 
The dynamical behavior of such process depends on the 
underlying structure of the system under study. Most of 
previous studies on coherent exciton dynamics are based 
on siinple structures, for example, the line 0, Q, cy- 
cle 0, Q , hypercube , Cayley tree [l^l , dendrimers [ll| , 
polymers [12| and other regular networks with simple 
topology. To the best of our knowledge, the dynamics 
of exciton on random network have not received much 
attention [isj . 

In this paper, we consider the coherent exciton trans- 
port on random networks of Erdos-Renyi (ER) . The co- 
herent exciton dynamics is modeled by continuous-time 
quantum walks (CTQWs), which is a quantum version of 
the classical random walk and widely studied by various 
researchers to describe the relaxation processes in com- 
plex systems El]- I^i the mathematical literature, 
the ER random network is defined as follows [l^ [13, [3 ■ 
Starting with N disconnected nodes, every pair of nodes 
is connected with probability p (0 < p < 1) and mul- 
tiple connections are prohibited. The ER random net- 
work is one of the oldest and best studied models of 
networks, and possesses the considerable advantage of 
being exactly solvable for many of its average proper- 
ties in the limit of large network size [19^. For instance, 
one interesting feature, which was demonstrated in their 
original papers, is that the model shows a phase tran- 
sition with increasing p at which a giant component 
forms (lol . [20| . An alternative and equivalent represen- 
tation of the ER random graph is to express the graph 
not in terms of p but in terms of the average degree k of 
the nodes, which is related to the connection probability 
p as: k = p{N — 1) « pN, where the last approximate 
equality is hold for large N. 

In the limit of large network size N, the degrees of ER 
random network follow a Poisson distribution peaked at 



the average degree k. In order to contrast the resem- 
blance and difference of the transport dynamics on net- 
works with the same average degree, we consider the co- 
herent exciton transport on a configuration model of ran- 
dom networks in which the degree k of each node equals 
to the average degree k {k £ Integers) of the Erdos- 
Renyi networks. The method for generating the graph is 
as follows [2l|: one assigns each node k (k £ Integers) 
"stubs" -ends of edges emerging from the nodes, and 
then one chooses pairs of these stubs uniformly at ran- 
dom and joins them together to make complete edges. 
When all stubs have been used up, the resulting graph 
is a random member of the ensemble of graphs with 
the equal degree [H, |2l|. The configuration model 
of random networks can also be implemented by using 
the edge-interchanging algorithm, which randomly inter- 
change two existing edges while keep the degree sequence 
unchanged [l^, |23| |. The configuration model is one of 
the most successful algorithms proposed for network for- 
mation, and has been extensively used as a null model 
in contraposition to real networks with the same degree 
distribution in biology, robustness, epidemics spreading 
and other dynamical processes taking place on complex 
networks |22l l24l . |25| . Here, we adopt this idea to com- 
pare the transport behavior on the two network models. 
As we will show, although the ensembles of ER model 
and configuration model have the same average number 
of connections, the transport dynamics on the two mod- 
els are different. 

The paper is structured as follows: In the next sec- 
tion, we briefly review the classical and quantum trans- 
port on networks presented in Refs. [1^, [23|. In Sec. Ill 
we study the time evolution of the ensemble averaged 
return probability on ER networks with different param- 
eters. Section IV presents the efficiencies of the classi- 
cal and quantum mechanical transport, and try to reveal 
how the model parameter affects the transport efficiency. 
In Sec. V, we consider the distribution of the long time 
averaged transition probabilities, and explore how the 
average return probability is related to network parame- 
ters. In Sec. VI, we consider the transport dynamics on 
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extremely connected networks. Conclusions and discus- 
sions are given in the last part, Sec. VII. 



II. TRANSPORT ON NETWORKS 

The coherent exciton dynamics on a connected net- 
work is modeled by the continuous-time quantum walks 
(CTQWs), which is obtained by replacing the Hamil- 
tonian of the system by the classical transfer matrix, 
H = -T The transfer matrix T relates to the 

Laplace matrix by T = —jA [10]. Here, for the sake 
of simplicity, we assume the transmission rate 7 for all 
connections equals to 1. The Laplace matrix A has non- 
diagonal elements Aij equal to —1 if nodes i and j are 
connected and otherwise. The diagonal elements An 
equal to the number of total links connected to node 
i, i.e., Aii equals to the degree of node i. The states 
\j > endowed with the nodes j of the network form 
a complete, ortho-normalised basis set, which span the 
whole accessible Hilbert space, i.e., J2k 1^ >^ ^1 ~ 
< k\j >— 6kj- The transport processes are governed by 
the master equation or Schrodinger equation [l^]- The 
classical and quantum mechanical transition probabili- 
ties to go from the state |j > at time to the state 
\k > at time t are given by Pkj{t) =< fc|e~*'*|j > and 
7Tk,At) = = I < fc|e-^*«|j > |2 0, respec- 

tively. Generally speaking, to calculate the transition 
probabilities, all the eigenvalues and eigenvectors of the 
transfer operator and Hamiltonian are required. We use 
En to represent the nth eigenvalue of H and denote 
the orthonormalized eigenstate of Hamiltonian by \qn >, 
such that J2n \ln X qn\ = 1- The classical and quan- 
tum transition probabilities between two nodes can be 
written as, 

n 

and 

-E„,e-*(^-^') (2) 
X < k\qn >< qn\j >< k\qi >< qi\j > . 

The above equations give the general expressions of the 
classical and quantum transition probabilities, which ex- 
plicitly depends on the eigenvalues and eigenvectors of 
the transfer matrix or Hamiltonian. A particular feature 
related to the transport is the return probability, which is 
the probability of finding the exciton at the initial node. 
The transition probability depends on the specific topol- 
ogy of the generated single network, therefore it is ap- 
propriate to consider its ensemble averages. 

III. AVERAGED RETURN PROBABILITIES 

The average of the classical and quantum return prob- 
abilities pjj (t) and TTjj (i) over all nodes of the network 
are 

Pit) = W E„ e'*^" < Qn\j ><j\qn> 



and 

= ^E„,e-*(^"-^') (4) 

X Ej < j\qn >< qn\j >< j\qi >< qi\j > ■ 

The classical p{t) is only dependent on the eigenvalues 
and decays monotonically from p(0) = 1 to the equipar- 
tition limt.^ooP(^) = 1/N. The quantum Tr(t) is depen- 
dent on the eigenvalues and eigenvectors, which is cum- 
bersome in the numerical calculations. The above equa- 
tions present the average of return probabilities over all 
nodes on a specific single network. In order to reduce 
the statistical fluctuation, we further average the return 
probabilities over distinct single networks, i.e., 

r— 1 

and 

where the index r denotes the rth generated ER net- 
work. Throughout this paper, we denote the average 
over network nodes by a bar (e.g., k, p{t), 7r(i), etc.), 
and the average over different realizations by a bracket 
(e.g., < p{t) >) while the actual values by undecorated 
characters. 

Figure [ija) shows the ensemble averaged return prob- 
abilities < p{t) > and < 7f(i) > on ER networks of size 
N — 100 with average degree k = 10, 20 and 30. For 
classical transport < p(t) > reaches the equipartition 
limt^oo p{t) = l/N very quickly. The curves at interme- 
diate times follow stretched exponential decay, which dif- 
fers from power law decay (t""-^) for the cycle graph [29j . 
The exponential decay of < p{t) > indicates that a classi- 
cal excitation will quickly spread the whole network and 
occupy each node with an uniform probability 1/A^ in a 
short time. It is evident that the excitation reaches the 
equipartition 1/N more quickly on networks with more 
connections (compare the curves in Fig.[lja)). For quan- 
tum transport < TT{t) > also decays quickly in the in- 
termediate times and then reach a final plateau. This 
plateau is larger than the equip-partitioned probability 
After a careful examination, we find such plateau 
corresponds to a constant value < •7r(i) 0.065 ± 0.01. 
Increasing the average degree k nearly does not change 
this value (compare the curves in Fig.[lja)). We note that 
here < n{t) > is smooth and does not display the oscilla- 
tory behavior, in contrast to the case for the cycle graph 
in which the return probability shows a quasi-periodic 
pattern [2§|. The non-periodic behavior of < 7f(t) > may 
be attributed to the large connectivity of the networks 
consider here. 

In Fig. [T{b), we show the same plot of < p{t) > and 
< 7f(i) > on a random graph model in which the degrees 
of each node are exactly equal to fc = 10, 20 and 30. The 
behavior of < p{t) > is almost the same as that on the ER 
networks. However, < Tt{t) > is quite different. < tt (t) > 
oscillates at intermediate times and also reaches a con- 
stant value. Such constant value (0.028 ± 0.003) is lower 
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FIG. 1: (Color online) Averaged return probabilities < p{t) > 
and < 7f(t) > for random networks generated by the ER 
model (a) and configuration model (b) with different values of 
degree. The networks are of size A'^ = 100 and the (average) 
degree of the networks are k — 10, fc = 20 and fc = 30 (see the 
different line styles) , respectively. All the curves are averaged 
over 100 independent realizations. 
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than that of the ER networks but larger than that of the 
cycle graph ((2iV — 2)/iV^ = 0.0198 for even-numbered 
networks and {2N - 1)/N'^ = 0.0199 for odd-numbered 
networks) [131 ■ The random networks generated by the 
ER model and configuration model have the same num- 
ber of total connections (on average), but the average 
return probabilities are quite different. Such difference 
may be caused by the different ensembles generated by 
the network models. 



IV. EFFICIENCY OF TRANSPORT ON 
INFINITE NETWORKS 

On finite networks, both the classical and quantum re- 
turn probabilities do not decay ad infinitum but reach a 
constant value at some time [231 ■ This value is related to 
the size of the networks. To reveal the decay behavior at 
large time scales, we consider < p{t) > and < 7f(t) > on 
infinite networks. In this case, the spectrum density can 
be regarded as a continuous distribution. Because the 
networks considered here are uncorrelated random net- 
works [3l[, the spectral density of the Laplacian Matrix 
converges to the semicircular distribution. 
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FIG. 2: (Color online) Efficiency of the classical and quantum 
transport. (a)Time evolution of p{t) and |Q(f)|^ for fc = 4. 
Both the classical and quantum transport at long time scales 
show the power law behavior. The exponent of local maxima 
of |Q(f)p is twice the exponent of p{t), implying the quantum 
transport is more efficient than the classical one. (b)Time 
evolution of p{t) and |Q(f)p for fc = 9. Here, p{t) drops 
below |ci(t)P and shows an exponential decay. Thus for k = 
9 the classical transport is more efficient than the quantum 
one. (c)|Q(t)p versus t for k — and k = 16. The curve 
depresses vertically on highly connected networks, suggesting 
that the quantum transport becomes more efficient when the 
connectivity of network increases. 



Where a = y^Np{l-p) and k = p{N - 1) for the ER 
networks [s^. This theorem is also known as Wigners 
law (33j , and it has extensive applications in statisti- 
cal physics, solid-state physics and complex quantum- 
mechanical systems [s^, llal ■ 

For sparse networks, i.e., p << 1, can be simplified 
as (7^ ~ k. Thus Eq. ([7]) can be written as. 



p(E) = < 2iTk ■ 

I 0, 



if \E-k\< 2Vk, 
Otherwise. 



(8) 



Therefore the spectral density is only a function of the 
average degree k on large sparsely connected networks. 
Using the above expression, we can calculate the return 
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probabilities of Eqs. 
as follows, 



([3]) and ^ in the continuum limit 



Pit) 
and 



p[E)dE 



2^ 



= \ J e-'*^ p{E)dEl 



e-'*g^4fc-(g-fc)^ 



dE, 
(9) 

(10) 



27rfe 



dE(^ 



Where the lower bound is obtained by using the Cauchy- 
Schwarz inequality and exact for regular networks [26l |. 
Analogous to the classical case, |Q:(i)p depends only on 
the eigenvalues of the Hamiltonian. Although is a 

lower bound and differs from the exact value 7f(i), |a (t)P 
quantitatively reproduce the overall behavior of tt (t) [29|. 
Therefore it is appropriate to use the decay of lafi)!^ to 
measure the efficiency of the quantum transport 2911. 

Fig. [2] shows the temporal behavior of and |Q!(t)p 
according to Eqs. ^ and pUj) for different values of av- 
erage degree k. For fc = 4 (Fig. [Ua)), p{t) scales as 
p{t) ~ t~^-^ and the local maxima of |a(t)p scales as 
|Q:(i)P ~ t^^ at long times. Such scaling argument can 
be understood by the spectral density. When fc = 4, the 
spectral density becomes as p{E) — ^/8E — E'^/8n. The 
long time behavior of p{t) and |Q;(t)p are mainly deter- 
mined by small E values, thus we can assume p{E) 
Ve [2^. Such scaling behavior leads to the power law 
behavior of the return probabilities, where the exponent 
for the quantum transport is twice the exponent of its 
classical counterpart (29|. 

The behavior of p{t) and |a(i)P alters accordingly 
when the average degree k increases. In Fig. ^h), we 
show p{t) and |Q:(i)p for k — 9. Now the classical p{t) 
does not show scaling but displays an exponential de- 
cay. We note that p{t) decays faster than |Q:(i)p, this 
indicates that the classical excitation spreads over the 
network faster than the quantum one. As a matter of 
fact, we find that for networks with fc > 4, p{t) always 
shows a faster decay at long times compared to |a(i)p. 
Such a behavior indicates that the classical transport is 
more efficient than the quantum transport on networks 
with large average degree. 

The behavior |Q!(i)P is qualitatively the same when the 
average degree k increases. To incarnate the difference 
of |q:(<)P, we plot |Q;(i)P as a function of t for fc = 16 
and fc = 64 in Fig. [UJc). As we can see, the local maxima 
of |a(i)P scales as |Q:(i)P ~ for both the values of 
k. However, the scaling depresses vertically for large k. 
This vertical depression of |a(t)p suggests that the quan- 
tum transport becomes more efficient when the network 
connectivity increases. 

We note that the increase of average degree greatly 
changes the efficiency of the classical transport. This is 
ascribed to large value of the eigenvalues E for network 
with large average degree. For large k, p{t) decays ex- 
ponentially, in contrast to the power law decay for small 
k. Quantum mechanically, the increase of average degree 
does not change the power law behavior of the local max- 
ima of |a(t)|2, but the first minimum is reached quickly 
and the local maxima of become lower on highly 

connected networks. 
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FIG. 3: (Coloronline) Ensemble averaged transition probabil- 
ity < Xk,3 > for random networks generated by the ER model 
(column (a)) and configuration model (column (b)) with dif- 
ferent values of degree 10, 20 and 30 (rows 1-3). The network 
size is TV = 100 and the average is over 100 realizations. The 
colormap is shown in the right hand of the plot. Dark regions 
denote large values of < Xk,j > and bright regions low values 
of < Xkj >. 



V. LONG TIME AVERAGES 

On finite networks, the transition probability con- 
verges to a certain value, this value is determined by the 
long time average. Classically, the long time averaged 
transition probability equals to the equal-partitioned 
probability 1/A^. However, the quantum mechanical 
transport does not lead to equipartition. Taking into 
account the ensemble average, we have, 

< Xk,j > (liniT-^oo ^ 7rk^j{t)dt) > 

(En I < %n >< qn\j >< k\qi >< qi\j > 



X lim7 



-it{E^ 



>< qn\j > 



=< {En,iS{E„-Ei)<k\q„ 
X < k\qi >< qi\j >) > . 

(11) 

Where d{E„ - Ei) ^ 1 for En = Ei and S{En - Ei) ^ 
else. Here, we numerically calculate the ensemble av- 
eraged transition probabilities according to the above 
equation. The results are show in Fig. [31 The dark 
regions denote large values of < Xk,j > and bright re- 
gions low values of < Xk,j >■ The dark blocks in the 
diagonal positions corresponds to the large return proba- 
bility. To view the quantitative behavior of < Xk,j >, we 
plot the distribution of < Xk.j > between all the pairs 
of two nodes in Fig. 21 The peaks in Fig. SKa) corre- 
spond to the non-diagonal transition probabilities while 



5 




V 



(b) 



k=10 
k=10 
k=20 
k=20 
k=30 
k=30 



<y > 

'•kj 




FIG. 5: Average return probability < x > on random net- 
works generated by the ER model (marked by k) and config- 
uration model (marked by k). (a)< x > versus k on networks 
of size A'^ = 100. (b)< x > versus network size A*' when the 
degree equals to 50. 



VI. TRANSPORT ON EXTREMELY 
CONNECTED NETWORKS 



FIG. 4: Distribution of ensemble averaged transition prob- 
ability < Xk,j > in the region 0.005 << Xfe.j >< 0.015 (a) 
and 0.02 << Xfcj >< 0.1 (b). The k with and without bar 
presents the random graphs generated by the ER model and 
configuration model, respectively. The peaks in (a) are from 
the non-diagonal transition probabilities and the peaks in (b) 
corresponds to the return probabilities. 



the peaks in Fig. IHa) correspond to the return prob- 
abilities [l^. The solid symbols denote the results for 
the ER networks, and the hollow symbols denote the nu- 
merical results for the random graphs generated by the 
configuration model where the degrees are exactly equal 
to the average degree of ER networks. It is observed that 
mean value of the non-diagonal transition probability for 
the ER model is smaller than that for the configuration 
model (see Fig. HJa)). A contrary conclusion for the re- 
turn probability is demonstrated in Fig. [IJb). We also 
find that the increase of degree does not change the cen- 
tral values of the peaks (compare the different symbols 
in the figure). This suggests that the quantum transition 
probabilities may do not change greatly when connectiv- 
ity of the networks increases [iJI • In order to reveal how 
the return probability is related to the average degree, 
we show both the node and ensemble averaged return 
probability (< x >=< J2jXj,]/N >) in Fig. [5); a). As 
we can see, the average return probability < x > is al- 
most a constant value in a wide range of average degree 
k but increases drastically when the network approaches 
to be fully connected 13]. Here, the difference for the 
two random graph models is also visible. The saturated 
value < X > of the ER networks is larger than that of the 
random graphs generated by the configuration model. In 
addition, we have also studied the relationship of the av- 
erage return probability < x > the network size N, 
which is shown in Fig.[5Ub). It is found that the averaged 
return probability < x > is inversely proportional to the 
network size N (see the asymptote indicated in the fig- 
ure) . This means that the return probability of small-size 
networks is larger than that of large-size networks. 



As we have shown, the average return probability 

< X > increases rapidly when the networks approach 
to be fully connected. To study this issue in detail, we 
consider the transport on extremely connected networks. 
In such case, the connection probability p approaches to 
1 and the degree of the network is of the order 0{N). 
We construct the extremely connected networks by ran- 
domly removing a certain number of edges on the fully 
connected network. The resulting network is equivalent 
to the ER random network, and the algebra of remov- 
ing edges is much computationally cheaper than the ER 
random graph algebra. 

Here, we consider the transport on networks of size 
N = 100, and assume the number of removed edges is 
TO. Fig. [S] shows the average return probability < x > as 
a function of the number of removed edges. We note 
that < X > decreases to ~ 0.1 when only 4% edges 
(200 edges) are removed. As more and more edges re- 
moved, < X > decreases slowly and tends to reach the 
saturated value ~ 0.65. We find that the < x > ver- 
sus m can be well described by an exponential decay 

< X >'~^ g-o.oi4r?x region to < 200. When only 
a few edges are removed, < x > is close to the return 
probability of the complete network (or fully connected 
network). For a fully connected network of size N, one 
eigenvalue of the Hamiltonian is and all the other values 
are equal to N, the long time averaged return probability 
is X = (N^ -2N + 2)/N^. This is a striking feature of 
CTQWs which differs from the classical counterpart. 



VII. CONCLUSIONS AND DISCUSSIONS 

In conclusion, we have studied the classical and quan- 
tum transport on Erdos-Renyi networks. We numeri- 
cally calculate the ensemble averaged transition prob- 
ability between two nodes of the networks. For finite 
networks, we find that the limiting transition probabil- 
ity is reached very quickly. For infinite networks whose 
spectral density follows the semicircle law, the efficien- 
cies of the classical and quantum mechanical transport 
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FIG. 6: (Color online) Average return probability < x > 
versus the number of removed edges m on networks of size 
A'' = 100. The straight line indicates the exponential decay 
for m < 200. 

are compared on networks of different average degree k. 
It is shown that the classical transport is more efficient 
than the quantum transport on networks with large con- 
nectivity. In the long time limiting, we consider the dis- 
tribution of the ensemble averaged transition probability, 
and show that the quantum transport exhibits a high re- 
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turn probability. Such high return probability almost do 
not change in a wide range of connection probability p 
but increases rapidly when the network approaches to 
be fully connected. For networks whose topology is not 
extremely connected, the return probability is inversely 
proportional to the network size N . In addition, we also 
compare the results with that on a random graph model 
in which the degree of each node equals to the average 
degree k of the Erdos-Renyi networks. 

We have shown that the transport dynamics on the ER 
networks is different from that on the equivalent network 
generated by the configuration model. The difference 
may be related to the different ensembles of the network 
model. Since the return probability reflects the symme- 
try of the network structure, the high return probability 
on ER networks may suggest a high symmetry of the 
topology [36j . 
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